© 2001 International Press 

Adv. Theor. Math. Phys. 5 (2001) 969-999 



Mathematical Analysis of 
the Photoelectric Effect 



Volker Bach, Frederic Klopp, and Heribert Zenk 

FB Mathematik; Johannes Gutenberg-Universitat; D-55099 Mainz; Germany; 
vbach @ mathematik. uni-mainz. de 

LAGA, Institut Galilee; Universite Paris-Nord; F-93430 Villetaneuse; France 
klopp @ math . univ -p ari s 1 3 . fr 

FB Mathematik; Johannes Gutenberg-Universitat; D-55099 Mainz; Germany; 
zenk @ mathematik. uni-mainz. de 

Abstract 

We study the photoelectric effect on the example of a simplified model of an atom 
with a single bound state, coupled to the quantized electromagnetic field. 

For this model, we show that Einstein's prediction for the photoelectric effect is 
qualitatively and quantitatively correct to leading order in the coupling parameter. More 
specifically, considering the ionization of the atom by an incident photon cloud consist- 
ing of N photons, we prove that the total ionized charge is additive in the N involved 
photons. Furthermore, if the photon cloud is approaching the atom from a large dis- 
tance, the kinetic energy of the ejected electron is shown to be given by the difference of 
the photon energy of each single photon in the photon cloud and the ionization energy. 
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I Introduction 

The photoelectric effect was discovered in increasingly precise experiments by Hertz [17, 
1887],Hallwachs [16, 1888], Lenard [18, 1902] 1 , andMillikan [21, 20, 1916]. 

It was observed that, when light is incident on a metal surface, electrons are ejected 
from the surface. The striking fact about this phenomenon is the seemingly odd dependence 
of the maximal kinetic energy T max of the ejected electrons on the frequency of the light 
and its independence of the light intensity. The latter contradicts the principles of classical 
physics, and in 1905, Einstein suggested an explanation of this phenomenon [9, 1905] 
which explicitly involves the quantum nature of electromagnetic radiation. He found that 

T max — hv — AE , (1.1) 

provided that the frequency v of the light times Planck's constant h, i.e., the photon energy 
hv is larger than the (material dependent) work function AE. Conversely, if 

hv < AE, (1.2) 

then no electrons leave the metal surface. Our ultimate goal is the derivation of Einstein's 
predictions, Eqs. (1. 1)— (1.2), from first principles of quantum mechanics and quantum field 
theory. 

In the present paper, we analyze a simplified model which is far from the appropriate 
model for a metal interacting with electromagnetic radiation and can, at best, be regarded 
as a caricature of a hydrogen atom interacting with the radiation field. Yet, it contains many 
of the mathematical difficulties we expect to encounter in the analysis of a more realistic 
model, and we prove Eqs. (1. 1)— (1.2) for this simplified model. Our emphasis lies in the 
following aspects: 

• Given the model as described in Sect. 1.1, below, our derivation is mathematically 
rigorous, and no unjustified approximations are used. To our knowledge, the present 
paper is the first to treat the photoelectric effect with mathematical rigor. 

We draw from many facts about nonrelativistic quantum electrodynamics which have 
been previously established in [2, 3, 4, 5]. 

As we propose to study the charge transported to infinity (see Sect. 1.2) which in- 
volves the asymptotics of the unitary time evolution operator e~ ltHt >, as t — > oo, 
our results can be viewed as part of scattering theory for models of nonrelativistic 
quantum electrodynamics. Results in this context, but on other aspects can be found 
in [1,6, 10, 12, 13,23]. 

• While our model for the particle system, i.e., the metal or atom, is a crude model in- 
volving only a single, spinless electron, the particle system is coupled to a quantized 
scalar field. (The difference between the quantized [vector] electromagnetic and a 
quantized scalar field is irrelevant, for the scope of this work. For certain other facts 
in nonrelativistic quantum electrodynamics, however, this difference is crucial, see, 
e.g., [4].) 

'in some physics textbooks, e.g., [15], Lenard is not mentioned, and it seems that in 1905, Einstein derived 
his famous theory from nothing but a Gedankenexperiment. 
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• Our model describes a single atom rather than a metal or even a gas interacting 
with the electromagnetic field. Thus our derivation of Eqs. (1.1)— (1.2) shows that the 
photoelectric effect is not a collective, statistical phenomenon, visible only if many 
particles or many photons are involved. Most texts on laser theory or quantum optics, 
e.g., [8, 14], immediately proceed to a statistical description of both the metal or gas 
of atoms, say, and the photon field, and the question, whether this is really necessary 
or merely a matter of mathematical convenience, is left open. 

• Within our framework, we prove Eqs. (1. 1)— (1.2) to be correct in leading nonvanish- 
ing order in the coupling parameter g which, in appropriate units, equals ^/2tt c?I 2 , 
where a m 1/137 is the fine structure constant. More precisely, given the interacting 
atom at rest plus an incident photon cloud consisting of N photons, we show that to 
leading order the contribution to the charge ejected from the atom is additive in each 
photon and independent of all other photons of the incident photon cloud. Moreover, 
we prove that this contribution is in accordance with (1. 1)— (1.2) in case that the pho- 
tons are in an incoming scattering state. In fact, the leading order contribution to the 
ejected charge resembles the first term in the Born series for the T-Matrix (see, e.g., 
[22]). 

• While it is customary to restrict the analysis of the photoelectric effect to a single 
photon scattering off the atom or metal, we point out that it is important to consider 
more than one photon, N > 2, because for a photon state consisting of a single 
photon only, total energy and energy of each single photon involved agree. Hence, if 
studying a single photon state, it is impossible to say whether the ejection of electrons 
is proportional to the total energy of the photon cloud or depends on the maximal 
energy of all photons in the photon cloud. 

Our paper is organized as follows. In Sect. 1. 1 below we introduce the mathematical model 
for the photoelectric effect. This includes a precise description of the atom, the photons, 
and their interaction in terms of a semibounded, selfadjoint generator of the dynamics, the 
Hamiltonian. 

In Sect. 1.2 we then describe our main results. In Subsects. 1.2. 1-1.2.2, we introduce 
the main quantities dealt with in this paper, the charge transported to infinity and photon 
clouds, and in Subsect. 1.2.3, we derive the asymptotics of the former to leading order in 
the coupling constant. A limit of monochromatic light is discussed in Subsect. 1.2.4, and in 
Subsect. 1.2.5, we compare our methods and results to those derived or used in other papers 
on scattering theory. 

The proof of the theorems in Sect. 1.2 are given in detail in Sects. II and III. 

Finally, our paper contains two appendices. In Appendix A we construct a Bogoli- 
ubov transformation that eliminates a single, arbitrary matrix element in the interaction, 
and in Appendix B we show that bound or negative energy states do not contribute to the 
transported charge. 

Acknowledgement: We thank S. De Bievre, J. Frohlich, Ch. Gerard, M. Griesemer, V. Ko- 
strykin, T. Paul, B. Schlein, R. Schrader, I. M. Sigal, H. Spohn, and S. Teufel for helpful 
discussions and remarks. 



972 



Bach, Klopp, and Zenk 



We gratefully acknowledge the support of our project by the TMR network ERBFM- 
RXCT960001 of the European Union and the Gerhard Hess-grant Ba 1477/3-1 of the ger- 
man science foundation DFG. 

1.1 Mathematical Model for the Photoelectric Effect 

Our goal is to analyze the photoelectric effect on the example of an atom coupled to the 
quantized radiation field. The mathematical model for our analysis of the photoelectric 
effect is a variant of the standard model of nonrelativistic quantum electrodynamics as 
given in [2] which we briefly recall and adapt to our problem at hand. 

1.1.1 Atom with a Single Bound State 

For the atom to be described we make the simplifying assumptions that it consists of a 
single, spinless electron bound in a potential well which admits exactly one bound state of 
energy e < 0. More specifically, we assume the Hamiltonian generating the dynamics of 
the electron to be given in diagonal form as 

H el = -A®e = (" A ° ), (1.3) 

acting on 

n el ■.= l 2 (r 3 ) e c = n ac {H el ) © n d (H el ) , (i.4) 

and H sc (H e i) = {0}. We further denote by 

ft := ( I ° ) and «■=-*-(»;) <W 

the projections onto the continuous subspace and the discrete subspace (of dimension one), 
respectively. The Hamiltonian H e i can be derived from a Schrodinger operator — A — V(x) 
with a short-range potential V. An appropriate choice of V guarantees that —A — V(x) 
has a single bound state only, and then — A — V (x) and H e i are unitarily equivalent, as can 
be seen by conjugating — A — V(x) with the wave operators. 

1.1.2 Photons 

We couple the atom described above to the quantized photon field which, for notational 
convenience, is assumed to be scalar. For our study of the photoelectric effect, the differ- 
ence is not relevant. The Hilbert space T, carrying the photon degrees of freedom is the 
bosonic Fock space T = JF b [L 2 (R 3 )] over the one-photon Hilbert space L 2 (R 3 ), i.e., 

oo 

T = ^ (n) , (1-6) 

where is the state space of all n-photon states, the n-photon sector. The vacuum 
sector, is one-dimensional and spanned by the normalized vacuum vector, f2, i.e., 
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:= CO. For n > 1, the n-photon sector is the subspace of _L 2 [(M 3 )™] containing 
all totally symmetric vectors. The Hamiltonian on T representing the energy of the free 
photon field is given by 

H f = J d 3 kLu(k)a*(k)a(k) , (1.7) 

where uj{k) := |fe| is the photon dispersion, and a*, a are the usual standard creation- and 
annihilation operators on T representing the canonical commutation relations, [a(k) , a(k')] 
= [a*(k),a*(k')} = 0, [a{k), a*(k')} = 5(k - k'), a(k)Q = 0, in the sense of operator- 
valued distributions. 



1.1.3 The Atom-Photon System 

The Hilbert space of states of the atom-photon system is the tensor product space 

H :=H el ®T , (1.8) 
and its dynamics is generated by the Hamiltonian 

H g := H + gW, (1.9) 

where 

H Q := H el ®l f + l el ®H f (1. 10) 

is the non-interacting Hamiltonian, < g <C 1 is a small coupling parameter, and 

W := J d 3 k\G{k)®a*{k) + G*(fc)®a(fc)} (1.11) 

is the interaction operator. Here, G G L 2 [M 3 ; B(H e i)) is a square-integrable function with 
values in the bounded operators on Tt e i, given by 

G(k) := ( 5,5 P t ] ) . (L12) 



Pl(k) 

for k e M 3 , a.e. For the formulation of assumptions about G it is convenient to introduce 

J Kn {k) := max ||<9£G(fc)|| + II (|x| 7 1) G{k)\\ + \\G(k) (\xp ® . (1.13) 

Note that, for J^ i7 (A;) < oo, we necessarily assume a decay of the coupling matrix G(k) 
as least as l^l -7 , as x — > oo. This assumption is not satisfied, e.g., if we consider elec- 
trons minimally coupled to the quantized radiation field. Namely, in the case of minimal 
coupling, the linear part of the interaction operator is of the form (1. 11), with 

Bm.c.(k) = )^ 1/2 gW-Vs, (1.14) 

where k = (k, a) € R 3 x Z2, k is a smooth function of rapid decay, as \k\ — > 00, serving 
as an ultraviolet cutoff, and e(k, ±1) _L k are two transversal, normalized polarization 
vectors. Besides the lack of decay of B m , c ,(k), as |a;| — » 00, B m , c ,(k) is not bounded. 
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This is, however, only a minor complication. The extension of our results to this case of 
main physical interest will be the subject of a forthcoming paper. 

Further note that 

J K ,-y(k) < J' K Jk) := max \\d%B(k)\\ + II |x| 7 5(fc) II + llB(jfe) |x| 7 II (1.15) 

\a\<K 

+ max ||9£pf(fc)|| + max ||<9j*p|(fc)|| + II |cc| 7 pf(fc)|| + |x| 7 II , 

and, conversely, J' K (k) < 3J^ i7 (fc). We shall make use of the following Hypothesis 
throughout the paper. 

Hypothesis 1. There is an integer K > 1 and a real number 7 > 3/2 such that G G 
C K [R 3 ; B(H e i)} is K times differentiable and satisfies the following estimate, 

1 + u(fc)~ 1 ) \J K ^(k)\ 2 d 3 k < 1. (1.16) 

We remark that, although the requirement (1.16) with fractional derivatives of G, for 
real K > 1, is presumably sufficient, we do not try to optimize our result in this respect 
and work with classical derivatives in this paper. 

Since, for a.e. k E M 3 , 

p T (fc):C -» L 2 (R 3 ) and Pl (k) : L 2 (R 3 ) -» C, (1.17) 

there exist p( - , k),rj(- k) e L 2 (R 3 ), such that 

[pi(k)z](x) = z V (x,k) and Pl (k)^> = - f p(x,k)rb(x) d 3 x , (1.18) 

for all z € C, ip € L 2 (K 3 ), and (x, k) <E M 3 x M 3 , a.e. Furthermore, in many applications 
B(k) G ^[L 2 (R 3 )] acts in the Schrodinger representation as a multiplication operator. That 
is, there is a function M e L 2 (R 3 x M 3 ) such that 

[B(k)i>](x) = M(x,k)tp(x), (1.19) 

for (x, k) G M 3 x M 3 , a.e. For instance, in case of the dipole approximation, 

M dip (x, k) = k(x/R) n{k) \k\ 1/2 e{k) ■ x , (1.20) 

using the same notation as for B m . c ., above, and additionally a spatial cutoff n{x/R) at 
length scale R 3> 1 which should be chosen large, compared to atomic length scales. We 
remark that, for a given function M, a physically natural choice for 77 and p is 

r](x,k) := p{x,k) := M(x,k) (p el (x) , (1.21) 

where ip e i e L 2 (R 3 ), \\(p e i\\ = 1, is the normalized wave function of the atomic bound 
state. Note that the special form (1.12) of G(k) implies that 

(P d ® 1) W (P d ® 1) = , (1.22) 
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since the lower right matrix element of G(k) vanishes. 

We remark that if the electron Hamiltonian is a Schrodinger operator, as described 
below Eq. (1.5), a single diagonal matrix element can always be "gauged away" by means 
of a Bogoliubov transformation. The construction of this Bogoliubov transformation is 
given in Appendix A. That is, for our model with a single atomic bound state, we do 
not loose any generality by assuming (1.22). As (1.22) is an important assumption for the 
present paper, we point out that, in case the atom has more than one bound state, Pd would 
be the orthogonal projection onto these, and (1.22) would imply, that the field does not 
couple different bound states. 

1.1.4 Selfadjointness, Semiboundedness, and Binding 

Next, we discuss selfadjointness, semiboundedness and existence of a ground state (bind- 
ing) of the Hamiltonian H g of the system. We basically invoke the theorems and methods 
for their proof established in [2, Sects. II and III]. 

Assuming Hypothesis 1, selfadjointness and semiboundedness, for any value of the 
coupling parameter g > 0, are a simple consequence of standard Kato perturbation theory; 
we appeal to [2, Cor. 1.7 and Lemma 1.5] which yield that H g is selfadjoint on its natural 
domain, dom[iJ g ] = dom[7Jo], and that the ground state energy is given by 

E Q := Ma[H g ] > e + O(g 2 ). (1.23) 

In the present paper we additionally assume that Eq is an eigenvalue of H g and that the 
corresponding eigenvectors have a large component in the vacuum sector . More pre- 
cisely, we require 

Hypothesis 2. The Hamiltonian H g possesses a (normalizable) ground state $ gs e H, 
||$ gs || = 1, i.e., E is an eigenvalue with corresponding eigenvector <& gs , 

H g $ gs = E $ gs . (1.24) 

Moreover, denoting Pq := = 1 <g) the ground state <I> gs obeys 

\\Pn$ S s\\ < Cg, (1.25) 

for some constant C < oo. 

The existence (1.24) of a ground state $ gs and the overlap bound (1.25) is proved in [2, 
Thm. 1.1] under the assumption of a somewhat stronger bound than Eq. (1.16) in Hypothe- 
sis 1, namely, 

J (l+uj(k)- 2 ^j\J 0n {k)\ 2 d 3 k < 1. (1.26) 

In this paper, we do not use this estimate but only its consequence in form of Hypothesis 2. 

We remark that the noninteracting Hamiltonian H has the unique (non-degenerate) 
ground state © f2 € H corresponding to its ground state energy eo, and we note in paren- 
theses that, by the results of [2], Eq. (1.25) holds for any ground state of H g . Consequently, 
the interacting ground state $ gs of H g is unique, provided g > is sufficiently small. 
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In Lemma III.2 we strengthen (1.25) and show that, for any a > 1 and any A < oo, 

|| {Hf + l) a P c $ gs || + || (ffj A) + l) a P£ $ gs || < O(g) . (1.27) 

Here, Hj denotes a free photon Hamiltonian acting only on photon states with energy 
less than A < oo, 

H f A) = I d 3 kuj A (k)a*(k)a(k) = [ d 3 kuj(k) a*(k)a(k) , (1.28) 

J J{uj(k)<A} 

where uj A (k) := u(k) l {ul ( k )<A}- 

1.2 Main Results and Discussion 
1.2.1 Charge transported to Infinity 

Having introduced the quantum mechanical framework in terms of Hilbert spaces and 
Hamiltonians and established some basic facts such as selfadjointness, semiboundedness, 
and binding, we proceed to defining the charge transported to infinity or simply the trans- 
ported charge. To this end, we introduce the projection F R onto the functions with support 
outside the ball of radius R > 0. More precisely, 

F R := ( *> J ) ® 1/ , d.29) 

where 1{\ X \> R ] := 1r3\b(o,r)[ x ]- Similarly, we introduce the projection Tj- onto the 
particle states with momentum in a measurable set T C K 3 , i.e., the functions whose 
Fourier transform is supported in T. That is, 

T T := ( \ er > jj ) ® 1, . (L30) 

where l{ pe r} acts as a Fourier multiplier with lr[p]- Given a state $ e the corre- 
sponding (least and most) transported charges with momentum in T C M 3 are defined to 
be 

liminf liminf II T r F R e~ itH <i * || 2 , (1.31) 
limsup limsup || T T F R e"^ 9 * || 2 , (1.32) 

R ^oo t^oo 

In case that T = R 3 , we write Q inf (*) := Q£ f (*) and Q su p(*) := Q™ p (tf ). Interpret- 
ing \1/ as an initial state for t = 0, the transported charge Q lnl (<&) measures the amount 
of its mass that is definitely transported away from the atom as time evolves, eventually. 
If T C R 3 is a small ball, then QrW additionally filters out the part of the state with 
momentum in T. 

While it would be desirable, of course, to define only one type of transported charge, 
namely, QrW := lim,R->oo lim t ^oo \\T r F R e~' ltH n *|| 2 , our methods described below 



Qt (*) ~ 
Q s r up (*) := 
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do not allow us to prove the existence of such a limit - not even for the restricted class 
of initial states of the form (1.35). Nevertheless, Q™ l (^f) and Q™ p (5') agree to lead- 
ing order in the coupling constant g, as we demonstrate below. In fact, our main ob- 
jective is the determination of the transported charges Qlf 1 and Q™ p (^) to leading 
order in g. Finally, we remark that it is important to observe the order of the two limits 
R — > oo and t — > oo, as interchanging these limits would, indeed, yield the trivial result 
lim^oo Hindoo ||Fr e - itH « *|| 2 = 0, for all 

As we prove in Appendix B, 

lim wv\F R \ m {H a )tr* H >*\ = 0, (1.33) 
lim sup||Ffll R -(iJ ff )e- itff »tf|| = 0, (1.34) 

R^oo t>0 

for all * e H, so the transported charge of all bound states W e Ranl pp (i? 3 ) and of all 
states W g Ranl R - (H g ) of negative total energy vanishes. 



1.2.2 Ground State and Photon Cloud 

Our choice for the initial state is of the form 

* := A(tJ)% s , (1.35) 

where 

A(0,f) = A(f) := l eI ® o*(/i)a*(/ 2 ) ■ ■ ■ a*(/ w ) , (1.36) 

and 

A(tJ) := e-" H * e lTHa A{f)e- lTHa e lTH <> . (1.37) 

The operator A(t, /) is called a photon cloud and represents photons with corresponding 
smooth orbitals / = / 2 , . . . , /at), /j g C^°(R 3 \ {0}), of compact support away from 
zero momentum. These N photons are prepared at time t = — r in such a way that they 
hit the atom at rest at time t = 0. (We thank S. De Bievre, M. Griesemer, H. Spohn and 
especially S. Teufel for clarifying this point to us.) From Eq. (II. 8) below we show that, 
in particular, A(oo, /) ( f > g s := lim T ^oo A(t, f)Q gs exists, thus representing an incoming 
photon scattering state. In contrast, choosing ^> according to (1.35) and with t = amounts 
to adding at time t = the photon cloud A(f) to the (interacting) atom at rest. 

We henceforth often leave out trivial tensor factors in our notation whenever it is clear 
from the context what is meant. For instance, we write A(f) = a*(fi) a*(f 2 ) ■ ■ ■ a* (/at), 
H = H ei + Hf, etc. 



1.2.3 The Photoelectric Effect 

Having defined the transported charges with momentum in T, we can now formulate a 
quantitative assertion about it. To this end, we fix p e M 3 and define two distributions 
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€ £>' (R 3 \ {0}) on smooth functions / € (M 3 \ {0}) of compact support away 
from zero by 

(I& , /) := 2tt J S(p 2 -E - w(k)) p(p, k) f(k) d 3 k , (1.38) 

{L °' f) - J p>-E -u(k) + iO (09) 

= -TT y 5(p 2 - S - p(p, fc) /(fc) d 3 fc (1.40) 

_ [ p(p,k)f(k)d 3 k 
J p2-E - u(k) ' 

P J denoting the Cauchy principal value and p(p, k) denoting the partial Fourier transform 
of p(x, k) (as a square-integrable function) with respect to the particle position variable x. 

Now we are in position to formulate our main result. 

Theorem LI. Assume Hypotheses 1-2, fix N, m£N with m < N, and let ipi , ip 2 , ■ ■ ■ <p m € 
Cq° (E 3 \{0}) be an orthonormal system, {<fii\<Pj} — Sij. LetA(r,ip) be the corresponding 
photon cloud, i.e., 

A(r,<p) = e~ iTH ° e lTHa a*((^i) ni a*(tp 2 ) n2 ■ ■ ■ a* (tp m ) n ™ e^ 11 " e lTH * , (1.41) 
where tij G N are such that m + n 2 + . . . + n m = N, and fix a measurable set T C R 3 . 

(i) If t > g~ x ' K then the transported charges with momentum in T of the initial state 
A(t, (p)& g s satisfy 

Q?(A(t^)% s ) + 0(g 2+ ») = Q™P{A(T,ip)% s ) + 0(g 2+ ») (1.42) 

m p 

= g 2 Q T r (v) := g 2 (n 1 \n 2 \---n m l) ^nj J ^f^p . 



j=i 



( ii) If t < g^ then the transported charge with momentum in T of the initial state 
A(t, <p)Q gs is given by 

Q?(A(t,v)% s ) + 0(g 2+ ») = Q^(A(t,v)%s) + 0(g 2+ ») (1.43) 

m „ 

= 9 2 Q T r(<P) ■■= 9 2 (n 1 \n 2 \---n m \)^2 n j / \(L P Q , ^)\ 2 d 3 p . 

3 = 1 Jr 

Here, and L V Q are the distributions defined in Eqs. (I.38)-(I.40), and p = p{K) := 
1 — K^ 1 > 0, where K is the degree of differentiability in Hypothesis 1 and 2. 



The asymptotic representation (1.42) (resp. (1.43)) is similar to the Born expansion for 
the T-matrix (see e.g. [22]) in Schrodinger scattering theory (resp. to the Born expansion 
for the wave operator). 



Mathematical Analysis of the Photoelectric Effect 



979 



Note that the coefficient Q T r (^p) of the leading order contribution to the photoelectric 
effect is additive (counting multiplicities) in the single photon states ipi, . . . , (p m . This 
implies that the contribution of each ipj is independent of the other single photon states ipi, 
i j, in the photon cloud, in contrast to the classical model for the photoelectric effect 
which would have predicted a dependence of Q^-{^p) only on the total energy of the photon 
cloud. 



1.2.4 The Limit of Monochromatic Light 

For large values of r, i.e., in Case (i) in Theorem 1.1 above, the transported charge directly 
yields Einstein's predictions (1.1) and (1.2). In Case (ii), however, there is an inconsistency 
between Eq. (1.43) and (1.1) and (1.2), since (L^ , ipj) may be non-vanishing, even if ipj is 
supported in a region where uj(k) < p 2 — E n . This observation may reflect that Case (ii) 
in Theorem 1. 1 is physically less relevant than Case (i), because it can hardly be realized 
experimentally: the photons would need to be present at the origin x = at time t — 
0, coming from nowhere. On the other hand, Case (ii) deals with a (perhaps naive, but 
reasonable) first proposal for a model of a photon cloud, and it is appropriate to discuss this 
case, as well. 

It turns out, that Einstein's predictions (1. 1) and (1.2) can be recovered even in Case (ii), 
provided the incident light is sufficiently monochromatic, i.e., sharply localized in momen- 
tum space. To define our notion of monochromatic light 

in precise terms, we restrict ourselves to considering a photon cloud consisting of single 
photon with wave function ipg G C£°(R 3 \ {0}), only. It is a trivial matter to extend the 
consideration below to a photon cloud of N orthonormal photon wave functions. 

For fixed uj > 0, we construct a wave function ip$ localized in energy about uj by 
choosing a smooth, L 2 (R)-normalized function of one, compactly supported variable, 
Xi,X2i • ■ • Xm € Cq°[(— 5 , 5)], and a smooth, normalized functions on the two-sphere, 
k e C°°[S 2 ] and setting 

p 4 (r,fi fc ) := uj- 1 5- 1 ' 2 x{^—p-) <n k ) , (1.44) 

using spherical coordinates (r, Ofe) = (|fc|, k/\k\). The limit of monochromatic light 
is then defined to be 5 — > 0. Note that <pg is asymptotically normalized in this limit, 

lim^o \\ips\\ = 1. 

Given the photon wave function ipg, a (more or less straightforward) computation in 
distribution theory yields, for any measurable set T C M 3 , that the first order term of the 
transported charge in Case (ii) in Theorem 1.1 is given by 

„ ( if E +uj<0, 

lim / \{L P Q , w)\ 2 d 3 p = \ 9 /o if £o+^ = 0, (1.45) 
6 ^° J r [ ®oI+ if E +uj>0, 
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where, 



9 := / 1 t (V e o+u<t) ^/E a +oj \G(y/E~ + uj a , uj)\ 2 d 2 a , 
Js 2 



0(p,r) 



• p(p, ra) n(a) d 2 a , 



with p(p, k) denoting the Fourier transform of x — > fc), and 



/o := 



x(a;)c?a; 



/OO /-C 
-oo J — ( 



y — x + iQ 
X(x)dx 



y — x + iO 



dy . 



dy . 



d-46) 
(1.47) 

(1.48) 
(1.49) 



We state Eq. (1.45) without proof but mention that its derivation requires a slightly stronger 
assumption than what is provided by Hypothesis 1-2. 

Einstein's predictions (1.1) and (1.2) are now manifest in the right hand side of (1.45). 
The coefficient lim^o Qri^-S^gs) may still vanish, as the factor 9o may turn out to be 
zero. This case is in accordance with (1.1) and (1.2) and merely reflects the physical fact 
that some transitions (here: from a bound into a scattering state) are forbidden in lowest 
order perturbation theory. 



1.2.5 Comparison to other Results and Methods in Scattering Theory 

As pointed out in the introduction, the results of the present paper may be considered scat- 
tering theoretic, as the transported charge derives from the asymptotic limit of the unitary 
evolution operator e~ %tH <> , as t — > oo. Most of the work in mathematical scattering the- 
ory in the past two decades or so is devoted to proving asymptotic completeness (AC), 
first for nonrelativistic A -body systems and, more recently and inspired by the former, for 
quantum field theoretic models like the one defined by the Hamiltonian H g , especially, in 
[6, 7, 10, 11, 12]. The mathematical methods developed and applied in these papers are, 
perhaps, more sophisticated than those used by us for the computation of the transported 
charge. A central role in all the papers mentioned above is played by asymptotic observ- 
ables, notably the asymptotic velocity operator. Loosely speaking, the positivity of the 
asymptotic velocity guarantees that any initial state eventually breaks up into parts which 
move independently, i.e., whose time evolution is free (noninteracting), and this is an im- 
portant ingredient to prove AC. It is therefore important to note that the transported charge 
measures any contribution that escapes to infinity (in particle configuration space), even 
those that come from particles with zero asymptotic velocity. Another important aspect to 
note is that we are interested in a quantitative estimate on the transported charge, asymp- 
totically in g. In this context, AC has to be seen as an existence result: every vector can 
be approximated by a polynomial in asymptotic creation operators, acting on the ground 
states of the system. It is not easy to turn the involved estimates leading to this statement 
into quantitative information. Also, it should be noted that so far AC can be proved only 
for massive quantum fields - an assumption that plays no role in our analysis, although we 
do make a confinement assumption on the interaction couplings. 
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To be more specific about the implications of AC on the value of the transported charge, 
we consider the results in [10], where AC is proved for Rayleigh scattering, i.e., for the 
energy range below the ionization threshold. In [10] an auxiliary Hilbert space H = W <g) 
T is introduced, and with the scattering identification operator I : H — > H given by 
I[<p <2> a* (fi) ■ ■ ■ a* (/„)£!] := a*(/i) • • • a*(f n )<p, the dynamics generated by H := H g ® 
ljr + In ® Hf on H can be compared to the dynamics on H generated by H g . The main 
result in [10] is, that 

0+ := s - lim e lrH * I e - lT " {l p JH g ) ® ljr) (1.50) 

T — >00 

exists and Ranf2 + D Ranl R - (H g ). In particular, for each ^ e Ranl R - (H g ), there exists 
V? € 7Y, such that 

limsup limsup ||Fr e ~* tHs *| = lim sup lim sup ||Fr e -4 *^ 3 fi + = 

R— >oo i^oc R^oo t^co 

= limsup limsup 7 e _i * H (g) 1^)^|| = 0, (1.51) 

R^oo t^oo 

and the last equation follows from the fact that e~ ttH leaves l pp (7i) invariant, and the 
application of I does not affect the electronic component, thus (B.l) implies (1.51). Hence, 
the proof of AC in [ 1 0] yields Q inf = Q sup = below the ionization threshold. While this 
result does not cover the range of positive energies, which is where the photoelectric effect 
takes place, the proof of AC for Compton scattering in the more recent paper [11] may 
show that Q mi = Q sup extends to any bounded energy range. It would then be interesting 
to see whether the methods used therein also yield a better quantitative estimate on the 
transported charge than the one derived in the present paper. 



II Dyson Series and Photoelectric Effect 

In this section we study the asymptotics, as t — > oo, for e~ ttHs A(t, f)Q gs , where H g — 
Hq + gW is the Hamiltonian generating the dynamics of the interacting system, which is 
applied to an initial state A(t, /) < E > gs for < r < g 11 and for t > g~~ . The initial state is 
composed of a photon cloud 

A = A(f) = a*(f 1 )a*(f 2 )---a*(f N ), with h,h,...,f N e C^(B A \ {0}) , (HI) 

of N photons and its time evolution 

A(t)<S> ss = A(r,/)$ gs := e^ 7 "^ e lrH ° A(f) e' irH ° e irIi ^> $ gs (0.2) 

applied to a ground state vector $ gs , 

H g $ gs = E $ gs , (H.3) 

representing the interacting atom at rest. Note that each photon orbital fj is assumed to 
be smooth and compactly supported away from zero momentum and in the ball of radius 
A < oo about the origin, 

supp/j C B\ \ {0} := {k e M 3 \ < u(k) < A} . (H.4) 
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We now give a list of assumptions used to derive Theorem II. 1, below. Given an operator 
X on H and t € K we denote its free Heisenberg time evolution by 

X t := e - ttHo X e ltH " . (11.5) 

Next, we recall from (1.27) the estimate 

||(tfj A) + l) a P c $ gs || + \\{Hf + l) a P^ $ gs | | < O(g) , (11.6) 

for any a > 1 and any A < oo, where is defined in (1.28). Eq. (II.6) is proved in 
Lemma III. 2, below. Furthermore, we recall from (1.22) that 

P d WP d = (P d ® 1) W (P d ® 1) = . (H.7) 

In Lemma III. 3, we demonstrate that, for any set, 

\\[W, A S ](^} A) + 1)- 1 -W 2 )|| < 0((l + \.s\)- K ) , (II.8) 

where if > 1 is the degree of differentiability in Hypothesis 1. Since $ gs e dom[(ii/ + 
l) 1+(w/2) ], for any N e N , the estimate (II.8) and the fact that H f and P d commute 
imply the norm-convergence of 

^ x {A{f),T) := lim * t (A(/),t) , (0.9) 

i — >oo 



where 



$ T (X,r) := ^ 



/ ds e^ -- 60 ) [W , X s ] P d $ gs if < t < 
Jo 

T 

J dse^ H °- E ^[W,X s }P d ^ es if r> 5 -* 



(11.10) 



-T 



provided if > 1. Furthermore, we show in Lemma III.4 that, for any r€l and s € Rq~, 

U^e-'^PctW, ^ s ](F/ + l)- 2 - (JV/2) || < G((l + |r|)- 3 / 2 ) . (11.11) 
Under these assumptions, we prove the following Theorem. 

Theorem II.l. Assume Hypotheses 1-2 and (II.l)-(II.ll). Then there exists a constant 
C < oo such that, for all t > g~ 1/>K , 

y-it(H g -E ) A( T )$ gs - A t 4> gs + ige-** 1 * - 130 ) $oo(A,t)\\ < C g 1+ ^ K) , (11.12) 

provided that either < r < g^ or r > anc/ where /i(if ) := 1 — K^ 1 > one/ 

if > 1 w the degree of differentiability in Hypothesis 1. 



Proof. First we compute the Heisenberg time evolution of the photon cloud A{t) on Ran 
[(H g -i)-^ N ], 

e - itH n A(j^ e itH si — e -i(t+r)H g e i(t+r)H ^ e -i(t+r)H e i(t+r)H g (H.13) 

S = t + T 

A t + 



e -isH g e isH Q e -isH & isH 9 

= A t - ig f +T dse- tsH * [W,A t _ s ] e isH ° 
Jo 

= A t - ig J ds e-^-^o [W , A s ] e l ^ H ^> , 
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using A = A(0) = A(0, /) and a change of integration variables, s i— ► t — s. Applying 
(11.13) to the ground state $ gs , we obtain 

-iff*(t,r) := e-^ H o- E ^ A{r)<S> gs - A t $ gs (11.14) 

= -igj dse- i ^- s ^ H "- E ^[W , A s }<P gs , 

and our goal is to show that 

||tf(t,r) - e-^-^^oo^W.T))!! < Ctf . (11.15) 



To this end, we first prove (11.15) for the special case r = 0. Given T > 1 and if > 1, 
and assuming that t > T, we observe that 



|$oo(A0) - #r(i4,0)|| < ©(T 1 -^) , 



(11.16) 



< ©(T 1 -^) . (11.17) 



Next, we use that H and A commute with Pd and that Pd W Pd = which imply 



Hence we have that 



\ e ~ UH <> <f> T (A,0) - e- ltHo <5>t(A,0)\ 



(11.18) 



(11.19) 



9 



f dr e-"' H <> W e ~ l{t - r)Ho $ T (A, 0) 
Jo 

f dr f ds e- tsEo e-" H ° W e -*(*-»-<0»o p c [^ , A s ] P d $ 
Jo Jo 



'0 Jo 
< fflK^ + l^+W 2 )*. 



f dr f ds\\We- l< - t - r - s)Ho P c [W,A s }(H f + iy 
Jo Jo 

and (II. 11) and a change of variable r t — r now imply that 

\\e- itH ° $ T {A,0) - e- itH ° $t(A,0)\\ < dgf ds 



■2— (JV/2) I 



o jo (l + | r -a|)3/2 
< C l C 2 gT, (11.20) 



where C\ < 00 and C 2 := /^(l + M) 3 ^ 2 dr < 00. Applying the estimates (11.16)- 



984 



Bach, Klopp, and Zenk 



(11.20) to the left side of (11.15), we hence obtain that 

||*(t,0) - e -it(H -E ) $^(^0)11 



(11.21) 



< 



*(t, 0) - / ds e-Ht-'HBg-Eo) [ W , A a ] $ g 
Jo 

f dse-^-^o-^ [W, A s ]$ gs - e- u ^- Eo ^T{A,0)\ 
Jo 



+ 

+ \\( e -it(H 9 -E ) _ e -i«(-ffo-Bo)-) $T ( A)0 )|| + ||$ T (A,0) - $oo(A,0)| 

/ dse is ^- E °^ [W, A s ]<S> gs - 9 T (A,0) + 0(T 1 - K + gT) 
Jo 



< 



Jo 



(H g -Eo) _ p is(H -E ) 



)[W, A s ]P d <P 



+ 



\ ds e is ^~ E ^ [W , A 8 ] P c $ gs + 0(T 1 - K + gT) , 
Jo 



using that P d + P c = 1. It remains to bound the two integrals in the last line of (11.21). The 
second of these is bounded by 



f T 

/ dse is ( H °- E °) [W, A S ]P C % 
Jo 

< ||(tf} A) + l) 1+ W 2 )p c $ gs || £ ds\\[W,A s ](H^ + 



(11.22) 



(AT/2) I 



< O(g) < O(gT) , 



due to (II. 6) and (II. 8), since K > 1 and T > 1. To estimate the first integral, we use again 
the DuHamel formula as in (11.19). This yields 



Jo 

/ e- ts£o / dr e irH ° W e ^ s - r)H ° P C [W , A s ] P d $ g 
Jo Jo 



(11.23) 



< g\\(H f + l) 2+ l N M$ 



gs 



/ ds / dr ||VKe i(s - r)Ho P c [VK, A s ](iJ/ + l)- 
Jo Jo 

f ds f 
JO Jo 



■2 — (JV/2) I 



(1 + | r |)3/2 



O(gT) 



changing again the integration variable r ^ s — r. Inserting (11.22) and (11.23) into (11.21), 
we arrive at 



|*(t,0) - e-^""-^ $oo(A,0)|| < 0(T 1 - K ' + 0T) 



(11.24) 



which yields (11.15) upon choosing T := gV( K ) 1 and setting := 1 — K 1 > 0, for 

the special case r = 0. 
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Note that if more generally < r < then $oo (A, t) = (A, 0) and hence 

||A(r)$ gs - A$ gs || = g\\ f dse- is ( H °- E ^[W,A_ s ]$ es \\ . (11.25) 

Jo 

Additionally using (II.8), we hence conclude || ^4(-r)<I>g S - A$ gs || < 0(g 1+fl ), proving the 
asymptotic expansion (11.12) in the case t e [0, g^]. 

In order to deal with the case r e [fl 1-1 ^, 00), we observe that an application of (11.25), 
for t > g~*, in combination with the decay property (II. 8) establishes 

\\e- ttH *A(t)<P gs - e- ItH *A(T)<t> gs \\ < 0{g 1+ ») , (11.26) 

and we may replace the time-evolved state e %tHa A(r)$ gs by e~ ltH <> A(i)$ gs in this case. 
Using r = t in Eq. (11.13), we furthermore obtain 

e -it(H g -E ) A ^ gs = At$gg _ ig J dse -i( t -s)H g [w^A^e^-^ 11 ^^ . (11.27) 
Using 

tf(i,t) = y dae _i(t - i)(ff »-' Bo) [W, A s ] $ gs , (11.28) 
computations running along the lines of (II. 16)— (II. 23) substituting 3>oo(A 00), $t(A 00), 

T 

00) for $00 (A 0), <&t(A 0), 0) and replacing integrals of the form J ds by those 



T 

of the form J ds proves (II. 12) for r > g~~<? , as well. □ 

-T 

Theorem II. 1 is the main tool for the analysis of the charges transported to infinity, 
Q 1 ^ 1 (A(r)<I> gs ) and Q s ^ p (A(T)$ gs ), because it allows us to take the limits t — > 00 and 
R — > 00, as shown in the following lemma. 

Lemma II.2. Assume Hypotheses 1-2, and let A = a*(/i) a*(/ 2 ) ■ ■ ■ 

a *{fN) with /1, /at g C^°(R 3 \ {0}). 77ze« //zere ex/ifi a constant C < 00, such that 

the charges transported to infinity obey 

Q™ p (A(T)<P gs )-Cg 2 +^ < g 2 \\T T <S>oc(A,T)\\ 2 < Q£ f (A(r)$ gs ) + Cg 2 +^ (11.29) 

/or < r < g M and t > g~i, where fi(K) := 1 — K^ 1 > and K > 1 z's f/ze degree of 
differentiability in Hypothesis 1. 

Proof: First we observe that, since F R = F R ® 1/ commutes with A t = l e z <8> A t and 
i?/ = l e ; ® if/, we have 

1 1 2V -Fr A t $ gs 1 1 < \\F R A t <S> gs \\ = \\A t F R <S> gs \\ (11.30) 
< \\A t (H f + iy N / 2 \\ ■ \\F R (H f + l) N ' 2 $ gs || . 

Applying Lemma III. 1 iV-times, we observe that ||A t {Hf + l) _Ar ' 2 || is bounded by C2 < 
00, uniformly in t, 



N 

J/2 



\A t (H f + l)- N / 2 \\ < JJ (fl- / + l)«- 1 )/ 2 o*(e- ito / j )(H/ + l) 



< C 2 . (11.31) 
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Moreover, (H f + l) N / 2 $ gs e H and F R — > strongly, as i? — > oo, hence 

lim lim 1 1 Tq- F R A t <& gs 1 1 = 0. (11.32) 

According to Theorem II. 1 we have 

e -it(H e -E ) A ^ = _ jge-it^o-Bo)^^^ + Rem (f) ; (11.33) 

where |jRem(£)|| < C\g 1+tl , for some constant C\ < oo independent of t e R+ and 
/i = 1 - l/if. Inserting Eqs. (11.33) and (11.32) into the definitions (1.3 1)-(1.32) of the 
transported charges, we find that 

Q^ up (A(r)$ gs ) - C l9 2+ » < g 2 Q T (A(r)% s ) < Q£ f (A(r)$ gs ) + C l9 2+ » , (11.34) 

where 

g r (A(r)$ gs ) := lim lim llrrF^*^-^)*^^^))!! 2 , (11.35) 

and estimate (11.29) would follow from Q r (^(r)$ gs ) = ||Tr$oo(^4,T)|| 2 . 

Next, we turn to TTF R e~ it( - H °- Eo ^ 00 (A, r). We note that F R = F R P c and hence 

p Re -it{Ho-E ) = e -it(H f -E Q ) F R P c e~ itHel = e~ u{H f- Eo) F R e- lt( -~ A) P c . 

(11.36) 

Furthermore, the absolute continuity of the spectrum of — A on i 2 (R 3 ) implies that 

(1 - F R )e-^- A ^ = l {w<K} e-^- A > - (11.37) 

strongly, as t — > oo, for any R < oo. Using (11.36), (11.37), and the fact that the operators 

T r = T T P C , e - lt( -- A \ and e -^(H f -E ) commute5 we thus obtain 

lim lim 1 1 TV F R e' lt{H °- Eo) ^(A^Ml = (11.38) 
= lim lim ||T T F fl e- it <- A >P c $ 00 (A,T)|| = \\T t $oo(A,t)\\ , 

R^oo t—>oo 

which implies the claim. □ 
Next, we recall from Eq. (II.9) the definition of 

POO 

/ dse ts{H °- Eo) [W, A s ]P d <5> gs if0<r<^ 

''' [ ° (H-39) 

/ dse^ H °~ E °'>[W,A s }P d ^ gs if T>g~ir 

J —OO 

which we complement by the definition of the vector 

pOO 

/ dse ls( - H °- Eo '> [W, A s }P d (0® SI) if0<r<<^ 
C(At) := { ° (11.40) 

/ dse ts{H °- Eo) [W, A s }P d (0Q) SI) ifr>g-T< 

J — oo 

whose existence is guaranteed thanks to Lemma III. 3. Lemma II. 3 below shows that these 
two vectors differ by at most O(g), i.e., a replacement of the ground state $ gs of H g by the 
ground state ® SI of H in the definition of ^^(A, t) introduces only small errors which 
are negligible, as we shall see. 
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Lemma II.3. Assume Hypotheses 1-2, and let A = a*(/i) <x*(/2) ■ ■ ■ o*(/jv) wif/i /i, 
G Co°(R 3 \ {0}). 77?en there exists a constant C < oo, rac/z f/zaf 



<M4 r) - C(Ar)|| < C.g. 



(11.41) 



Proof: We first note that due to Lemma III. 3, we have the estimate 

||$oo(A,r) - #£(4 r)|| < 
< 



(11.42) 



(M A) + l) 1+ f (P d $ gs - (0 8fi)) ■ / d*|| [W, A s ](H { f A) + 

Jr 

(H { f A) + i) 1+ %(p d <t> gs - (oefi)) 



for some constant Ci < 00. Additionally using Lemma III. 2 and the decomposition 

P d $ gs ~ (0 © fi) = [PdPn^gB - (0 © 0)] + P d P f f $ gs , we observe that 



(^} A) +l) 1+JV / 2 (P d $ gs - (0©Q)) < O(g) 



□ (11.43) 



The final ingredient for the proof of Theorem 1. 1 is Lemma II.4, below. Note that due 
to the definitions (11.40) we may restrict ourselves to the case r = and r = 00. 

Lemma II.4. Assume Hypotheses 1-2, fix N,m £ N with m < N, and let ip\ , (f2 , • • ■ <£m € 
C^°(R 3 \ {0}) be an orthonormal system, (ipi\<Pj) = 6%,j. Let A = a* (ipi) ni a* ((p 2 ) n2 ■ ■ ■ 
a*((p m ) nm , where nj £ N are such that ni + n 2 + . . . + n m = N, and fix a measurable 
set T C M 3 . Then 

\\T t ^(A,t)\\ 2 = Q r T (<p) = (11.44) 

m p 

= (m!n 2 !---n m !) Vnj / |(L? , ^)| 2 d 3 p , 

wifft LP = LlJorr £ [O,^], and = L^Jorr £ [g- 1/K , oo), wnere L£ anc/L^ are 
rfe/ined in (I.38)-(I.40). 

Proof: We first write = a*(G) + a(G) and use the canonical commutation relations to 
obtain that 

[W, A.] = [a*(G) + a(G) , a^e"^)"^^^)™ 2 • • • a* (e~ is " <p m ) n ">] 
= [a(G) , a*(e-^^ i ra*(e-^ lf2 r^--a*(e-^ lfm ) n -] 

= J2 n A ( G \ e ^^) ® Ha^e-^^r*- 5 ^ I , (11.45) 
j=i ^ i=i ^ 

on dom[ii^ +Ar / 2 ], where (G|e~' sa Vj) G B{H e i) is a bounded operator acting on the 
electron variables given by 

(G|e- <Sfc V j ) := J d 3 ke- lsu ^ ip 3 {k)G*{k) . (11.46) 
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Moreover, for any s > 0, 

m 

e u,{ - H °- Eo ^(G\e- Uu ip j )® J|a*(e- i ' a V j ) ni -*-' (0©O) 

i=i 

m 

= e^^-^^Gje-^^)® Jla^e-^^)"*-^ (0 0) 

i=l 
m 

= e ia (- A - B °)(P c GP < j|e- isa '^)® JJo*^)"'-^ (0©O). 

= (fefflOj^Ha'y"'^!!, (11.47) 



i=i 



where the Fourier transform of tpj e L 2 (M 3 ) with respect to the particle coordinate x is 
given by 

4»(P) = / d3k e ts( > p2 - E °-^ k » p(p, k) ip 3 (k) . (11.48) 

Here, we used the fact that P d G(k)*P d = and hence G{k)*P d = P c G(k)*P d , for all 
k e R 3 . Moreover, P c G(k)*P d (i>®z) = z Pl (k)* 00, and this yields (11.48). Eqs. (11.45)- 
(11.48) imply that 



(11.49) 



Passing to the momentum representation for the particle variable and using that 

< m m \ m 

l[a*( Vi r- 5 ^n Y[a*( Vi r- 5 ^n) = Sj,f Hint- 6^)1, (11.50) 



we hence obtain 



\Tt*°(A,t) 



m „oo rn 

j=i J - T »=i 
ni!n 2 ! •••n m ! J d 3 p J 



dsipj,s(p) 



(11.51) 



It remains to evaluate the integral J^° t ds tpj tS (p) in (II. 5 1), which is an exercise in distribu- 
tion theory. To this end, we note that the K > 1-fold partial differentiability of p(x, k) with 
respect to k in Hypothesis 1 implies that p(p, ■ )ipj is two times partially differentiable with 
respect to k, and that its derivatives of order < 2 have bounded support away from k = 0. 
Since the phase in ipj, s (p) is non-stationary, away from k = 0, two times integration by 
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parts back and forth yields 

/oo poo r 

dsi) jtS {p) = J ds J d 3 ke ls ( p2 - E °-^ k » p(p,k)cp 3 (k) 

= f^r Jd s ke^ 2 - E °-^(iV.0[p(p,k)^(k) 
= lim J d 3 kj^ rf se is(p 2 -So-^(fe)+te) pfak)<pj(k) 

. ie i(p 2 -E -u(k)) fifa fc ) 



(1 / d 3 fc^l 



■ Pj(fc) . 



o 2 -E - u(k) + is 
Thus in the case r = 0: 

/•OO />1 />OC 

/ ds^ Jy8 (p) = I dsip j)S (p) + / dsi> j<s (p) 
Jo Jo Jl 



(11.52) 



(11.53) 



s^oj p 2 - E a - Lu(k) + is J \ 



For the case r = oo, we get 



lim / dsi> js (p) = lim / ds [ d 3 k e ^(p 2 -Eo p ^ k ) = (n.54) 

^°°J-t t^°°J- t J 



-t 

r e it(p 2 -E -uj(k)) _ e it{p 2 -E -Lu{k)) 

lim / d 3 k — 2 — -j-^ p(p,k)ipj{k) 



i(p 2 - E - Lu{k)) 



= 2n [ d 2 fl(p 2 - E ) 2 p(p;p 2 - E ,n)^(p 2 - E ,fl). □ 

Proof of Theorem LI: The proof of Theorem 1. 1 consists of the following chain of esti- 
mates, 

Q^(A(t)% s ) + (D(g 2+ n = Q s ^(A(r)% s ) + 0{g 2+ ») 

= g 2 \\T T ^oo(A,T)\\ 2 (11.55) 
= g 2 \\T T ^(A,r)\\ 2 + 0(g 2 +n (11.56) 
= Q T T (<p) + 0(g 2+ n , (11.57) 

where Lemma II. 2, II. 3, and II. 4 justify Eq. (11.55), (11.56), and (11.57), respectively. □ 



III Technical Estimates 



In this section we derive the estimates (II. 6), (II. 8), and (II. 1 1), which is the basic input for 
the asymptotics of the time evolution asserted in Theorem II. 1 . 
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Before we turn to these bounds, we derive a preparatory lemma. 



Lemma DLL Suppose that a > 1, A < oo, and let G e L 2 [R 3 ; B(H)] with C := 
/ d 3 k(l +Lu(k)- 1 ) \\G{k)\\ 2 < oo. Denote a*{G) := / d 3 kG(k) ® a*(k) and a(G) := 
/ d 3 kG{k)* ®a(k). Then 

||( J ff} A) + l) a a(G)( J ff / + l)- a - 1 / 2 || < VC , (IH.1) 

||(ffJ A) + l) Q a*(G)(i7 / + l)- Q - 1/2 || < VC(2 + 4A) Q , (IH.2) 

||(i/} A) + l) a a(G A )( J tfj A) + l)" Q " 1/2 || < VC, (IH.3) 

||(fr} A) + l) a a*(G A )(H} A) + l)~ a ~ 1/2 || < VC(2 + 4A) Q , (ffl.4) 

where G\(k) := l{ w (fc)<A} G(fc), anJ l{ w (fc)<A} denotes the characteristic function of 

{k e R 3 | w(fc) < A}. 



Proof. First, we introduce a more compact notation, 



Gi .- G, Hi .= H f , ui .- w, 

r< n u u( A ) (III.5) 

Cr 2 := Ga, rl 2 := iiy , ^2 := W A , 

and observe that Eqs. (HE.l)— (III.4) are equivalent to 

||(ffj A) + l) Q a(G,)(if, +1)-"- 1 / 2 !! < VG, (m.6) 
||(F} A) + l) a a*(G j )(i?j + l)" a " 1/2 || < VG(2 + 4A) a , (III.7) 

with j = 1,2. We apply the operator on the left side of (III. 6) to a normalized vector tp e W 
and obtain the desired estimate by means of the pull-through formula, the Cauchy-Schwarz 
inequality, and HjV + 1 < ifj + 1 + Wj, 



(ijj A) + ira(G,)(^ + l) 



< 



< 



d 3 k \\Gj(k)\\ 



H 



-a-l/2 
(A) , i 







(III.8) 



a(fc) (flj- + l)- 1/2 ^ 



(/^ 



iiG,(fcr^y /2 
(fc) ) 



H; 



Hi + 1 



1/2 



< 



^(fc)|W | 1/2 
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for any tp € H, \\il>\\ = 1. To derive (III. 7), we use the canonical commutation relations, 



(Hf> + l) Q a*(Gj) (Hj + I)-"" 1 / 2 V|| (HI.9) 
= J M % (, | WGtf) « ( Hj + 1)— N*) « + 1)- «•(*) 

= J d 3 kd 3 k(a{k){H j + l)- 1 / 2 tlj 



(ff} A) + l+W A (fc)+WA(fe)) 



a(fc) ( J ff J + l)- 1 / 2 ^ 

.(^■ + l+^(fc))(^ + l + ^(fc))J 



< 



Hf ] + 1 + 2A 




" a (Jd 3 k\\G 3 (k)\\ 





+ 



(Hj + I)" 2 "" 1 (ff^ + 1 + C A (fc)) 2a || ( / d 3 k ||G,(fc)|| 2 ) 



< C(l + 2A) 2a (l+ / d 6 kwj(k) a(fc) (flj- + 1) _1/2 V 



< C(l + 2A) 2a 1 + 



if,- 



< C(2 + 4A) 



2a 



Note that in the last step we used i?j = J d 3 k ujj(k) a* (k)a(k). 



□ 



Now we are in position to establish Eq. (II.6) which shows that the interacting atom- 
photon ground state $ gs is well-localized in energy, even for the noninteracting Hamilto- 
nian. 



Lemma III.2. For any a > 1, A < oo, and sufficiently small g > 0, we have 



Hf ] + l) a P c $ g 



4 A) +irP f f<i> gs || < 0{g). 



(III. 10) 



Proof. The proof is similar to the one for [2, Thm. X]. We first note that the asserted bound 
(III. 10) is implied by (1.25) and the following estimate, 



{Hf> + 1) a 1 {H a >E /3} *gs|| < 0(g) 



(in. ii) 



because on Ran l{H <E a /3}, we have P. c = an d #y < l e ol — l-Eol/3. 

Let x € Cq°(/; [0, 1]) be a real-valued, smooth function, compactly supported in J := 
(|£7 , §£o) and such that X (E ) = 1. Clearly, 1 {j?0 >.e /3} x(#o) - 0, and X (#g) 3> gs = 
$ gs . Hence we observe that 



(h^ + iri {Ha > Ea m $ g 



(in. 12) 



{ h >e q/3} (H^ + l) a [ X (H g ) - x(H )] (H g -E Q + l)-"- 1 ^ . 
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Next, we represent x{H g ) and x(i?o) by the functional calculus based on almost analytic 
extensions of smooth functions of compact support [19]. Let x € Cq°(A; C) be an almost 
analytic extension of x, supported in a small complex neighborhood A C C of I and such 
that &x(z) = 0((lmz) 2 ). By means of the measure d/i(z) := (2Tri)~ 1 dx{z) dz A dz and 
the second resolvent equation, we have the following identity 

= -g J 'd»(z){(H -z)- 1 W(H g -z)- 1 } . 

Inserting this identity into (III. 12) and applying Lemma III. 1 to W = a*(G) + a(G), we 
obtain the norm estimate 



l(4 A) +iri { H > Eo /3 } $ g! 

< 9 



< 



(A) t 

J d»(z){(H - z)- 1 (Hf ] + l) a W 

■(Hg-Eo + iy-^Hg-z)- 1 }^ 
\d»(z)\{\\(H - z)- 1 !) ||(i/j A) + l) a W{H f + I)-"" 1 ! 

■\\(H f + ir +i (h 9 -e + I)— 1 !) \\(H g z)- 1 !)} 

< 0(g), (HI.14) 
which finishes the proof. □ 

Our next goal is the derivation of (II. 8). Actually, we prove a somewhat stronger es- 
timate which is also used in the proof of Lemma III. 4, below. Recall that A = A(f) = 
a*{h)a*{h) ■ ■ ■ a*(f N ) is a photon cloud of N photons /i, f 2 , . . . , f N G Of (B A \{0}), 
that X t := e~ ltHo Xe ltH ° denotes the free time evolution of an observable X, and that 
K > 1 is the degree of differentiability in Hypothesis 1 . Denoting 

(1 + N) 3 / 2 :=( (1 + I ^ I)3/2 !)«!, (BUS) 
we prove the following lemma. 

Lemma III.3. For any N g N and K > 1, there exists a constant C < oo such that 

||(l + |a:|) 3 / 2 [W,i4 8 ](H} A) + l)- 1 -W 2 )|| <C(l + \s\y K , (IH.16) 
for all set 

Proof. First observe that on dom(H^ /2 ), we have e~ isHo a* (f)e tsH ° = a*(e- lsu f) and 
hence 

A s = 8 '(e-' s 7i)«*(e- su / 2 )... fl '( e - su / ff ). (IH.17) 
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Writing W = a*(G) + a(G) and using the canonical commutation relations, we thus have 

[W,A S ] = [a%G)+a{G),a%e- is "h)a*{e- is "h)---a*{e- is "f N )] 

= [fl(G),a*(e-»/ 1 )a'(e- isu / 2 )...fl*(e-»/ w )] (111.18) 

N , N n 

where (G\e~ lsw fj) € B(H e i) is a bounded operator acting on the electron variables given 
by 

<G|e-" w /j> := J d 3 ke- isaj W fj{k)G*{k) . (01.19) 

Note that fe fj(k)G*(k) is K > I times continuously differentiable, and, thanks to the 
support properties of fj, there exists rj > such that fjG* G C*[£?a \ _B rj ; B(W e ;)] nas 
compact support away from zero. Thus, thanks to Hypothesis 1, K times (B(H e i) -valued) 
integration by parts yields the standard estimate for oscillatory integrals, 

(l + |a;|) 3/2 (G|e- isfc 7 j ) (IH.20) 

< \s\~ K J d 3 k ||(1 + \x\f/ 2 (V fc • {k/\k\)) K fj(k) G*(k)\ 

< o((i + \ s \)- K ) , (in.2i) 

provided \s\ > 1. For \s\ < 1, Estimate (111.20) is trivial and hence holds for all s £ 
R. Inserting (III. 20) into (III. 18) and undoing the free time evolution (which is possible 



because H commutes with H^), we obtain 



(1 + M) 3/2 [W,A S ] (Hf> + 1) 



(A) 



-1-(JV/2)| 



(in.22) 



N 



N 



< 0((1 + \s\)- k )Y: (Ua^e-^fi) (i/f'+l)" 1 ^) 

N , N s 

7=1 v .=i, ' 



-(JV/2) 



To estimate a product of TV — 1 creation operators 
/i 5 /2, • • ■ , 

/jv € C$°(B A \ {0}), we use Estimate (III.4) N — 1 times and derive 



o*(/i)o*(/ 2 )"-o*(/ JV _i)with 



JV — _L \ 

n «*(/«)) ( J ff} A) +i)- (jv - i)/ 

V-l , 

n ( (i?j A) + i)^- 1 )/^*^) (hj a) + 1)-/ 2 
• n IN A) + i) (n ^ i)/2 «*(/n) + 1)- /2 



(111.23) 



71=1 

N-l 



71=1 



< 0((1 + A)^ >) . 
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Inserting this bound into (III. 22) yields (III. 16), since A < oo is bounded. □ 
We finally turn to proving (II. 1 1). Note that the following Lemma uses Hf rather than 

<>. 

Lemma III.4. There exists a constant C < oo such that, for any r <G K and s G Rq , 

|| We- irHo P c [W, A s ](H f + l)- 2 - (N ^\\ < C ((1 + \r\)- 3 / 2 ) . (111.24) 

Proof. The proof is based on the decay of the Schwartz kernel 

e-^- A Hx,y) = (47r«r)- 3 /2 

cxp[i(x — y) 2 (4r) _1 ] of the propagator e~ lr (~ A ^ of the free particle. More precisely, 
defining C\ :— (Att)^ 1 j(l + |x|) 27 d 3 x < oo, we observe that 



1 



(in.25) 



I" 



xd A y 



< 



( 47rr )3/2 

\w\\ u\\ 

( 4?rr )3/2 yj (1 + | X |)27 



d 3 X 



< Ci 



-3/2 



for any tp,ip e L 1 C\ L 2 (R 3 ), since 7 > 3/2. This estimate yields 

||(l+|.T|)-7 e -^op c(1 + | a .|)-7|| 



(111.26) 



(l + |a;|)-Te- ir (- A )(l + |a;|)-T 




-irH i 



< dr- 3 / 2 . 



Inserting (III. 26) into (III. 24) and using that Hf and H commute, we obtain 

II W e- irHo P c [W , A s ] (H f + l)- 2 -( JV / 2 ) || 



(111.27) 



< \\w(i + \x\y ®(H f + i) 



-1/21 



|(1 + \x\)-'< e- lrIia P C {1 



||(i + \ x \y ® (Hf + 1) 1 ' 2 [w , a s ] (^ + i)- 2 -w 2 )|| 

< C ir -^ 2 \\W(l + \x\y® (Hf + l)- l ' 2 \\ 

||(1 + \x\y ® (ff, + l) 1 / 2 [W , A a ] (Hf + l)- 2 -( Ar / 2 )|| . 

Next, we use (III. 18) and the pull-through formula to commute (Hf + l) 1 / 2 through 
[W , A s ], observing that supp/j C £?a(0). This yields 



(1 + \x\r ® (if, + l) 1 / 2 [1+ , A s ] (i?/ + l)- 2 - W 2 ) 



(111.28) 



< 



(1 + \x\r [W , 4 a ] (H f + l)- 2 -(^/ 2 ) + 1 + (N - 1)A) 



1/2 



< (1 + A^A) 1 / 2 (l + |a;|) 7 [W, ^(Fj + l)- 3 / 2 -^/ 2 ) 



< C 2 , 



for some constant C2 < 00, uniformly in s e R. Here we additionally inserted (III. 16) to 
derive the last inequality. Finally, there exists a constant C 3 < 00, such that 



W(l + \x\y® (Hf + l)- 1 ' 2 



< C* 



(111.29) 



as follows from writing W (1 + |x|) 7 = a*[G (1 + |x|) 7 ] + a[G (1 + \x\)i] and applying 
Lemma III. 1 . □ 
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A Elimination of Matrix Elements by a Bogoliubov Trans- 
formation 

Consider H et = — A— V(x), a Schrodinger operator on R 3 having a single, simple negative 
eigenvalue, e < 0, associated to a normalized vector ip . This Hamiltonian represents an 
atom which we couple to the photonic field described in Subsect. 1.1.2 with an interaction 
described in Subsect. 1.1.3. More precisely, we assume the interaction W to be of the form 
(1. 11), with G G L 2 [R 3 ; B(H e i)] being a square-integrable function with values in the 
bounded operators on H e i, satisfying 

J d 3 k{l + tu(k)- 2 } \\G{k)\\ 2 < 1. (A.l) 

Note that this assumption is stronger than (1.16). We now show 

Theorem A.l. Assume (A.l). For g > small enough, there exists a unitary transform 
U g € B(H e i <8> T) real analytic in g such that 

H g := U g Hg U* g = H el ® 1 f + 1 el ® H f + gW (A.2) 

where 

• H e i = H e i + g 2 AV, and AV is a bounded self-adjoint operator; 

• ifW acts as a multiplication operator in the electron variable, AV is a potential, 
i.e., a multiplication operator; 

• H e i has a single, simple negative eigenvalue, e g < 0, associated to the one-dimensional 
eigenprojector H g ; 

• we have 

(n s ® i/) w (n g ® i f ) = o. (A.3) 

All the quantities introduced above are real-analytic in g. 

Theorem A. 1 proves that, starting from a Schrodinger operator on M? with a single, 
simple negative bound state, using a unitary transform, one can always pass to an interac- 
tion fulfilling Assumption (1.12). 

Proof. Pick h e L 2 (R 3 ), and consider the Bogoliubov transform 

U{h) = l^e'sKW+'Wl . (A.4) 

A standard computation gives 



U(h)H g U(h)* = H el ® 1/ + \ el ®H f + gW, 



(A.5) 
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where 

H el = H el + g 2 AV(h), (A.6) 
AV(h) := j d 3 k{h{k)G{k)+h(k)G*{k)} , (A.7) 

W = W + l®{a*(w/i) + a(w/i)}. (A.8) 

We remark that AV : L 2 (M. 3 ) — > B(H e i) is a bounded linear operator with norm bounded 
by one, thanks to (A.l). Hence, i/ e ; has a single, simple eigenvalue e(ft) = e + 0(g 2 ) in 
a vicinity of e , and we denote the corresponding normalized eigenvector by <p g (h). 

To prove Theorem A. 1, we only need to show that, for g > sufficiently small, we can 
construct h g £ L 2 (R 3 ) such that 

-(<P g (h g ),G{k)<p g {h g )) = h g (k), where G(k) := Lu(k)- 1 G(k) . (A.9) 

By Assumption (A.l), G is a square-integrable function with values in the bounded opera- 
tors on Hei ■ Using standard perturbation theory, we construct ip g (ft) simply by normalizing 
the vector 

g {h) := ^-j^ {z-H el -g 2 AV(h)Y\ dz, (A.10) 

where 70 is a small circle of center e and radius | e /2 1 . The function ip g (h) is real analytic 
in g and continuous in h, for h in the unit ball in i 2 (R 3 ); real analytic perturbation theory 
and the bound ||AU(/i)|| < \\h\\ immediately give 

\\<p g (h) - <p g (.h')\\ < Cg 2 \\h-h'\\, (A.ll) 

for some constant C < 00 and all ft, ft' £ L 2 (R 3 ) with ||ft||, ||ft'|| < 1. Set 

ft - -(<Po\G(k)<po) ■ (A.12) 
For g fixed and small and ft £ L 2 (R 3 ), with ||ft|| < 1, we define 

T g (h) = -(<p g (h + h)\G(-)cp g (h + h)) - ho. (A.13) 

Note that To has a trivial fixed point, h = 0. For g > sufficiently small, T g maps the unit 
ball of L 2 (R 3 ) into itself. Moreover, one computes 

||T s (ft)-T s (ft')|| < 4||G(fc)||||^(ft + ft)-^(fto + ft')ll 

< Cg 2 \\h-h'\\, (A.14) 

for some C < 00. Hence, T g is contracting for g small enough. Therefore, the fixed point 
equation h = T g (h) has a unique solution in the unit ball of L 2 (R 3 ). This fixed point is 
the desired function solving (A.9). We may construct this solution as the norm limit of the 
sequence T™(0). Each of these terms being real analytic in g and the convergence being 
uniform in g sufficiently small (the rate of convergence is given by g 2 ), the limit is real 
analytic in g. □ 
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B Transported Charge for Bound and Negative Energy 
States 

Theorem B.l. Assume Hypothesis 1. Then 

lim sup \\F R e- ltH * lpp(H g ) * II = (B.l) 
lim sup \\F R e- itH ° l,^]^) *|| = (B.2) 

R^oo t > 

for any \t € Ji, where \ pp (H g ) and lj^oo ojC-Hg) are the spectral projections of H g onto 
its point spectrum and onto ] — oo, 0], respectively. 

Proof. To prove equality (B.l), we first assume that ^ is an eigenvector of H g with corre- 
sponding eigenvalue E £ R. Since F R — > strongly, as i? — > oo, we then have 

lim sup||F K e^* H ' * II = lim ||Fr*|| = 0. (B.3) 

ft— >oc £>q R^oo 

This, of course, generalizes to any finite linear combination of eigenvectors of H g . Since 
the norm closure of this set of vectors is the pure point subspace of H g , Eq. (B.3) also 
generalizes to vectors in Ran l pp (H g ), which yields (B.l). As a consequence of (B.l) we 
may assume — l c (-ff s )\E' in the proof of (B.2). These vectors out of the continuous 
subspace have the property that 

lim l (E -6,E+6)(H g )* = 0, (B.4) 

for any E e R. As H g is bounded from below and in view of (B.4), it hence (choosing 
E := 0) suffices to prove that 



R 



lim \\F RX [H g }\\ = 0, (B.5) 



for any smooth function X <= Co° [( — 00 > 0)]' compactly supported on the negative half-axis 
and away from zero. The basic idea of our proof of (B.5) is essentially the same as [2, 
Thm. II. 1] or Lemma III. 2 and uses a Combes-Thomas or Agmon Estimate. For H = H g 
or H = H , we use a representation 

based on an almost analytic extension X <= Cg°(A4;C) of x, whose compact support 
M C C can be chosen to include only zeC with Rcz < —5, for some 5 > 0. Moreover, 
we can choose X as to obey Bx(z) = ©((Imz) 2 ). The measure in (B.6) is then defined as 
d/i(z) := (2-7ri) _1 Bx(z) dz A dz. By means of (B.6), the second resolvent equation, and 
the fact that F R x[Ho] = F R P c x[Hq] = 0, since H > on Ran P c , we thus have 

F RX [H g ] = F R (x[H g ] xlHo}) 

= - 9 [ d»(z)F R ^ c — w ir—- (R7) 

Jm -A + Hf-z H g - z 
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Next, we pick 1 < A < oo and introduce fx : M 3 — > R+ by fx (x) := y/l + (x/X) 2 , and 
we denote f\{x)^ 1 := 1 / fx{x). We note the following properties of fx, 

y < fx(x) < l + \x\ 3 ' 2 and fx(x) _1 \Vfx{x)\ < A -1 , (B.8) 
which imply 

\\F R f xix^W < CXR- 1 , (B.9) 

Re{/ A (x)(-A)/ A (x)- 1 } = (B.10) 
-A + 2Rc{/ a (.t)- 1 V/a(xHV} > -A" 2 , 

\\fx(x)W(H g +i)- 1 \\ < C, (B.ll) 

for some constant C < oo. Here, (B.9) is trivial, (B. 10) is meant in the quadratic form sense 
on smooth functions of compact support, and (B.ll) uses Hypothesis 1 and Lemma III. 1 . 
Choosing A > 1j\fl and inserting the bounds (B.9)-(B.ll) into (B.7), we arrive at (B.5). 

□ 
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